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1 Introduction 

By a resistor network, or simply network, we mean here a grapli Q = 
(y, E) witii vertices V and (undirected) edges E, and c a positive real- valued 
function on the edges. The value Ce is the conductance of the edge e. 

The Laplacian on a network (^,c) is the linear operator A: MY — t- M^ 
defined by 

{M){v) = Y^c^Af{v)-f{v')). (1) 

Here the sum is over neighbors v' oi v. 

The Laplacian is the most basic operator on a network, and has countless 
uses in all areas of mathematics. We study here some specific properties of 
the Laplacian on networks embedded on surfaces, and in particular on the 
simplest surfaces: the plane, the annulus, the torus. 

Our goals will be to focus on a few topological, probabilistic and combi- 
natorial applications. Specifically, our three main interconnected goals are 
to discuss: 

1. The discrete EIT (Electrical Impedance Tomography) problem: re- 
constructing the network from boundary measurements. For planar 
networks, the classification due to Yves Colin de Verdiere of Dirichlet- 
to-Neumann operators [Gj [7] , 

2. The connections with the random spanning tree model developed by 
Curtis, Ingerman and Morrow [8] and Wilson and myself [191 [2Q], and 

3. The characteristic polynomial of A on an annulus and torus based on 
joint work with Okounkov and Sheffield [16], Okounkov [15], and in 
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The main new results in the current paper are Lemma [4. Ill and Theorems 
EBEIlandO 

2 Background 

2.1 Electrical Impedance Tomography 

The classical Calderon problem is to determine the electrical conductivity 
of a region in M^ or M'^ by making boundary measurements of the type: fix 
the voltage along the boundary and measure the resulting current flow out 
of the region. 

This problem has concrete applications in medical imaging, nondestruc- 
tive materials testing, etc. In these settings it is called Electrical Impedance 
Tomography (EIT). It is widely used because it is easy and cheap to imple- 
ment: basically it is very easy to apply a low voltage to pretty much any 
conductive material and measure the resulting current flow. 

In this paper we study the discrete version of EIT. Given a resistor net- 
work (a finite graph in which every edge is a resistor), understand its internal 
structure from boundary measurements. There is a very nice solution for 
planar networks. For networks on more general surfaces we are beginning to 
have a better understanding but much more remains to be done. 

For a network ^ on a surface with nontrivial topology there is closely 
related question, even when there is no boundary, which is to understand the 
determinant of the "bundle Laplacian" as a function on the moduli space of 
fiat connections on bundles on Q. 

We shall see how these questions are related below. 

2.2 Spanning trees 

One of the combinatorial tools involved in the study of discrete EIT is the 
model of (random) spanning trees, and related objects called cycle-rooted 
spanning forests on networks. Kirchhoff [22] was the first to see the con- 
nection between spanning trees and the Laplacian on a network: he showed 
that the Laplacian determinant in fact counts spanning trees. 

The random spanning tree measure (or UST, for Uniform Spanning Tree) 
has for the past 20 years been a remarkably successful and rich area of study 
in probability theory. The random spanning tree as a probability model 



goes back to at least work of Temperley [21], who gave a bijection between 
spanning trees on the square grid and domino tihngs. Aldous and Broder 
[H 12] in the 1980s showed a connection between spanning trees and random 
walk, giving an algorithm for sampling a UST based on the simple random 
walk on the network. Pemantle |28] showed that the unique path between 
two points in a UST has the same distribution as the loop-erased random 
walk (LERW) between those two pointql|. Wilson [32] extended this to give a 
quick and elegant method of sampling a uniform spanning tree in any graph. 
Burton and Pemantle [4j proved that the edges of the uniform spanning tree 
form a "determinantal process", that is there is a matrix T indexed by the 
edges such that the probability that some set of edges {ei, . . . , 6^} is in the 
UST is det(T(ej, Gj)i<i,j<k)- The matrix T is built from the Green's function 
of A. 

The UST on the grid Z^ received particular attention due to the confor- 
mal invariance properties of its scaling limit (limit when the mesh size of the 
grid goes to zero). In particular in [13] we showed that the "winding field" 
of the UST on Z^ (which measure the winding of the branches around faces) 
converges in the scaling limit to the Gaussian free field, a conformally invari- 
ant random function on M^. In [25] Lawler, Schramm and Werner showed 
the convergence of the branches of the UST on Z^ in the scaling limit to the 
conformally invariant process SLE2. 

As mentioned above, Temperley [3T] gave a bijection between spanning 
trees and the dimer model on the square grid; this was later generalized to 
arbitrary graphs in p/Tj and [18], pointing to a very close relation between the 
UST and the dimer model. Many of the results discussed below are inspired 
by the corresponding results in the dimer model. 

2.3 Characteristic polynomials 

To a network on an annulus we can associate a polynomial P{z) which is 
the determinant of the line-bundle Laplacian A for a flat line bundle with 
monodromy 2; G C* around a generator of the homology of the annulus. 
One can think of this as the determinant of the Laplacian acting on the 
space of 2;-quasiperiodic functions: "locally defined" functions whose values 



^The loop-erased random walk between a and b is defined as follows: talk a simple 
random walk from a stopped when it reaches b, and then erase from the trace all loops in 
chronological order. What remains is a simple path from a to b. 



are multiplied by z under analytic continuation counterclockwise around the 
annulus. 

The polynomial P is real and symmetric (that is, Pi^z) = P{l/z)). There 
is a nice characterization of those P which occur: they are those with positive 
constant coefficient, and whose roots are real, positive and distinct except 
for a double root at z = 1. See Theorem 16.11 below. 

For a network on the torus one can associate a two-variable polynomial 
P{zi,Z2) in a similar way: it is the determinant of the Laplacian for a flat 
line bundle with monodromy zi, Z2 around generators of tti of the torus. This 
polynomial plays an important role in understanding the structure and large- 
scale properties of the random spanning tree model on a periodic planar graph 
[U [TB]. Like in the case of the annulus there is a characterization of those 
P[zi,Z2) which occur: the curve {P{zi,Z2) = 0} is a Harnack curve which 
is symmetric {P{zi, Z2) = P{zi^, 2;^^)), satisfies P(l, 1) = and has positive 
constant coefficient. See the definition of Harnack curve in section [731 One 
can show that in fact any Harnack curve satisfying the above properties arises 
from the characteristic polynomial of some toral network. 

A more recent development in networks on the torus is the understanding 
of the integrable dynamical system underlying them. Here the dynamical 
system is not related to the flow of current in the network, but is a dynamical 
system on the space of conductances assigned to the edges in the network. 
It is the restriction of an integrable Hamiltonian dynamical system on a 
larger space (the space of line bundles on bipartite networks on the torus) 
to a Lagrangian subvariety. See [12]. We will not discuss this here, mostly 
because we don't yet have a very good understanding of this system except 
as an invariant subvariety of a larger system. 

For networks on other surfaces S one can also define a characteristic 
polynomial P = det A which is a regular function on the representation 
variety of vri(S) into SL2(C). Here much less is known about the structure 
of P. See however Theorem 18.11 below. 

Finally, one would like to consider a network on a general surface with 
boundary, and understand properties of the boundary measurement map 
(for the appropriate bundle Laplacian) and to what extent it determines the 
network. We are still far from this goal but the cases of simple topology or 
small boundary hint that there is an interesting answer. 



3 Networks 

By a network we mean a connected finite graph Q = (V, E) witli a positive 
real- valued function c on tlie edges called the conductance function, and a 
usually nonempty set of vertices A/" called nodes or terminals or boundary 
vertices. These are the simplest type of electrical networks, with resistors 
only. 

The nodes are indexed from 1 to n. Let / be the set of non-node vertices; 
elements of / are referred to as internal vertices. 

3.1 Laplacian 

The Laplacian on a network Q is the operator A : MY — ?■ MY defined for 
/ G M^ by ([1]). Since MY has a natural basis, consisting of functions S^ 
which are 1 at f and zero elsewhere, we often write A as a matrix A = 
{Ayyi)yyifzv where A^yi = —c^v' if f,f' are adjacent and the diagonal entries 
are A^j^, = Ylv'r^v ^™' ■ ^^^ ^^^ example ([8]) which is the Laplacian of the graph 
of Figure [H 

Associated to the Laplacian is the quadratic form Q on MY defined as 

Q{f) = (/,A/) = J2c.Afiv)-fiv')r (2) 

which is called the Dirchlet energy of /. Here (,) is the natural inner 
product in which the basis {6y}^^v is orthonormal. 

From the expression ([2]) we see that A is positive semidefinite with kernel 
consisting of the constant functions (we assume unless stated otherwise that 
Q is connected). 

Another useful operator is the incidence matrix d, defined as follows. 
Fix an (arbitrary) orientation for each edge e, so that e+ and e_ are its head 
and tail. Define d: R^ — t- R^ by the formula 

ci/(e) = /(e+)-/(e_). 

Then it is easy to verify that 

A = d*Cd, 

where C is the diagonal matrix of conductances, and d* is the transpose of d 
(its adjoint for the inner product (, ) above): for an element w G R^ we have 



d*w{v) = y ^w(e) 



where the sum is over the edges e containing endpoint v, and directed towards 
V. We can call df the gradient of the function / and d*w the divergence 
of the flow w. 

Note that we can write the Dirichlet energy as 

Qif) = {df,Cdf) 

where here (,) is the natural inner product on M.^. 

3.2 The Laplacian with boundary 

li B G V is a subset of vertices and I = V \ B, we can define an operator 
Ab : M.^ ^ M.^ called the Dirichlet Laplacian with boundary B by the 
same formula 

ABfiv) = Y,c..'ifiv)~fiv')), 

where here f G / and the sum is over all neighbors of f in \^ (not just those 
in/). 

As a matrix A^ is just a submatrix of the full Laplacian A: the sub matrix 
indexed by /. 

3.3 The Dirichlet problem 

A function / G M^ is said to be harmonic at v if A/(f) = 0. By ([1]) 
this simply means that f{v) is a weighted average of the values of / at the 
neighboring vertices, where the weights are proportional to the conductances. 

A function is harmonic if it is harmonic at all vertices. On a finite 
connected graph a harmonic function must be constant, since the kernel of 
A consists of constant functions only. 

Suppose i? is a nonempty set boundary vertices and u is a function on 
B. A function fonV which is harmonic on I = V\B and agrees with u on 
B is said to be harmonic with Dirichlet boundary conditions on B. 
Such functions satisfy a maximum principle: their maximum (and mini- 
mum) values occur on B (because of the above observation about weighted 
averages) . 

Given this fact it is not hard to see that, if B is nonempty, then A^ is 
invertible: any function / in the kernel of A^ can be extended to a function 
on V taking values zero on i?; it is thus a harmonic function with Dirichlet 
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boundary conditions at B. By the maximum principle it must take its max- 
imum and minimum on B, where it has value 0. Thus / must be identically 
zero. 

The Dirichlet problem is the problem of finding, for a given function u 
on B, a. function / on ^ which is equal to u on i? and harmonic on I = V\B: 
find / satisfying 

Af{v) = OioTveV\B 

f{v) = u(v) for V E B. 

The function / is called the harmonic extension of u. For finite networks 
the solution to the Dirichlet problem is straightforward using invertibility of 
A^. Alternatively, the solution / is the function with boundary values u 
and which minimizes the Dirichlet energy Q{f). Convexity of Q implies that 
there is a unique solution. 

A function on Q which is harmonic on / is a solution to the Dirichlet 
problem with boundary B. Such a function is also called an equilibrium 
potential. 

3.4 The Dirichlet problem and the response matrix 

For a network we will take A/", the set of nodes, to play the role of the 
boundary vertices (which we called B above). Given a function u on A/", 
let fu be the harmonic extension of u. Then A/u(t>), restricted to A/", is a 
(generally nonzero) function of f G A/". The map A : -u H- A/u|_v is a linear 
map A : M-^^ — )■ IS/'^^ called the Dirichlet-to-Neumann map of the network. 
For certain conventional reasons we often prefer to work with L = —A and 
refer to L as the response matrix. See the next section for an example. 

We can describe A as a matrix as follows. Suppose we index the vertices 
so that vertices in A/" come first. Then we can write A in block form as 

where A has rows and columns indexed by the nodes A" and C has rows and 

columns indexed by the internal vertices; C is just the matrix of the Dirichlet 

Laplacian A^. 

The Dirichlet problem 02]) can be written as follows: find // on the internal 

vertices / so that 

'A B\ (u\ ^(g- 

B' C [fj U 



for a function g which is the desired function of u. We can spht this into two 
hnear equations 

Au + Bfi = g (4) 

B'u + Cfi = 0. (5) 

This second equation can be solved as 

// = -C-'B'u. (6) 

(Here C is invertible, as discussed above). 

One can then plug in this value of // in the first equation to get 

Au - BC-^B*u = g, 

that is, Am = g where A is the matrix A = A — BC^^B^ . Thus 

L = -A + BC-^BK (7) 

We see from this expression that L is symmetric, since A and C are 
symmetric. In fact A is positive semidefinite (and L is negative semi definite) , 
with kernel consisting solely of the constant functions. This follows from the 
fact that 

(m,Am) = (/,A/)>0 

where / is the harmonic extension of u (one can check this using (jH])). 



3.4.1 Example 

Let Q be the network of Figure [1] with three nodes {fi,f2,f3}, one internal 
vertex ^4 and three edges with conductances ci, C2, C3. We have 

ci -ci 

C2 -C2 

C3 -C3 

-Ci -C2 -C3 C1 + C2 + C3 



A 



(8) 



Using formula (jTj) we have 
ci 

- I C2 
C3 




Ci + C2 + C3 ' 



-C2 -C3) 



-C1C2-C1C3 

C1+C2+C3 
C2C1 

C1+C2+C3 

C3C1 
C1+C2+C3 



C1C2 



C1C3 



C1+C2+C3 
-C2C1-C2C3 

C1+C2+C3 

C3C2 
C1+C2+C3 



C1+C2+C3 

C2C3 

C1+C2+C3 

-C3C1-C3C2 

C1+C2+C3 



(9) 
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• vi 



V3 

Figure 1: The 'Y' network. 

3.5 Electrical interpretation 

Associated to a function / (called potential or voltage) on a network is a 
flow u on the edges (called current) defined by 

a;(e)=Crf/(e) = Ce(/(e+)-/(e_)). 

Here C is the diagonal matrix of edge conductances and df{e) = /(e_|_) — /(e_) 
is the gradient of /. 

Ohm's law states that Cdf{e) is the flow across edge e when the vertices 
are held at potentials /(e+) and /(e_), and the resistance on the edge is 
1/ce, the reciprocal of the conductance. 

Kirchhoff's network equations state that when the boundary nodes are 
held at fixed potentials, the internal vertices will attain potentials / in such 
a way that the flow of current induced by Ohm's law in the network has 
the property that the current entering an internal vertex equals the current 
exiting that vertex. In other words the divergence d*u of the current cu is 
zero at internal vertices. Thus d*Cd{f) = at every internal vertex, so / is 
the harmonic extension of u. 

This allows us to give Lij the following electrical interpretation: Hold 
all nodes except i at potential 0, and node i at potential 1. Then for the 
resulting current flow, Lij is the current exiting the network at node i. This 
current is positive: the potentials at internal vertices are in [0, 1] (in fact in 
(0, 1)) by the maximum principle for harmonic functions, so the current at 
edges adjacent to j is directed towards j. Similarly one shows La < 0. 

11 



4 Circular planar networks 

By a circular planar network (CPN) we mean a connected finite graph 
Q = (V, E) embedded in the plane, with a positive real-valued function c on 
the edges called the conductance function, and a nonempty set of vertices 
A/" on the outer face called nodes or terminals. (There may be other vertices 
on the outer face as well). The reason it is a circular planar network is that 
the network can be drawn inside the disk with the nodes on the boundary 
circle. 

The nodes are ordered in counterclockwise order from 1 to n. Let / be the 
set of non-node vertices; elements of I are referred to as internal vertices. 
The non-exterior faces of Q are referred to as interior faces; we consider 
Q to also have n exterior faces, with the ith exterior face comprising the 
region between nodes i and i -|- 1, and having boundary consisting of the 
edges along the "outer face" of Q between node i and i + 1. 

4.1 Conjugate harmonic function 

If ^ is a circular planar network, a harmonic function f on Q with Dirichlet 
boundary conditions on JV has a conjugate harmonic function g defined 
(up to an additive constant) on the faces of Q by the condition Cdf = dg, 
where d is the difference operator on the dual graph Q': if e = vv' is an 
edge of Q and Fi, F2 are the adjacent faces, so that Fi is left of e when e is 
traversed from f to f', then 

Ce{f{v')-f{v))=g{F,)-g{F2). (10) 

The harmonicity of / implies that these equations are consistent around 
a vertex; by summing along paths one sees that the value of g is defined 
consistently on all faces, once one choose a value of g at any starting face. 
Thus these equations define g uniquely up to a global additive constant. 
The function g is harmonic for the operator d*C^^d. That is, one puts 
conductances on the dual edges which are the inverses of the conductances 
on the primal edges. 

The pair (/, g) is sometimes written a.s f + ig and is called a discrete 
analytic function since the equations ( !T0|) play the role of the discrete 
Cauchy-Riemann equation. 

If / is harmonic with Dirichlet boundary conditions then g will be har- 
monic on Q' with Dirichlet boundary conditions, that is, g is harmonic on 
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interior vertices of Q' (that is, on non-exterior faces of Q). 

4.2 Characterization of response matrices 

The response matrix is a fundamental matrix associated to a network. A nat- 
ural question is: what properties of the network are reflected in its response 
matrix? 

The basic properties of L shown in the previous section are that 

1. L is real and symmetric, 

2. L is negative semi definite, 

3. The kernel of L consists of the constant functions, and 

4. Lij > if z 7^ j, and La < 0. 

(It is not hard to see that property 2 is implied by the others). 

Any matrix satisfying the above properties is in fact the response matrix 
of some (not necessarily planar) network: it suffices to take a complete graph 
on n vertices, declare that each vertex is a node, and put conductance Lij 
on the edge ij connecting vertices i and j. 

The question is much more interesting in the case of circular planar net- 
works. In this case there is a beautiful characterization of response matrices 
in terms of minors, due to Colin de Verdiere [B] . 

Given two disjoint subsets A, B of the nodes of a circular planar network, 
we say that A and B are noninterlaced if they are contained in disjoint 
intervals of nodes with respect to the circular order. That is, if a < 6 < c < o? 
are nodes then it is never the case that a,c & A and b,d & B. 

Theorem 4.1 ([6]). L is the response matrix of a connected circular planar 
network if and only if L is symmetric, with kernel consisting of the constant 
functions, and for any disjoint noninterlaced subsets of nodes A and B with 
\A\ = \B\, we have detL;^ > 0, where L^ is the suhmatrix of L whose rows 
index A and columns index B . 

There are two parts of this theorem, the necessity of the conditions and 
the existence of a network for any L satisfying the conditions. Both have 
interesting proofs which we will discuss. 

We define a network to be well-connected if all noninterlaced minors 
det L^ are strictly positive. In [6j Colin de Verdiere shows that the space of 
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well-connected circular planar networks on n nodes is liomeomorphic to an 
open ball of dimension I r, ) ■ More about this in section 14.4.21 

4.3 Groves 

As one can see from (jTj), the response matrix entries are rational functions 
of the conductances. Surprisingly, the off-diagonal entries are subtraction 
free rational functions of the conductances, that is, ratios of polynomials 
with positive coefficients; in fact the coefficients are all 1 or (see for example 
Q). This is a strong indication that these polynomials are counting some 
type of combinatorial objects. 

Indeed, that is one of the ideas in the Curtis-Ingerman-Morrow proof (see 
[8]) of necessity in Theorem 14.11 above, which we discuss in this section. 

4.3.1 Spanning trees 

A spanning tree of a graph is a subset of edges which is a tree, that is, con- 
tains no cycles, and is spanning, that is, connects all vertices. Equivalently 
one can say that a spanning tree is a maximal set of edges (with respect to 
inclusion) that has no cycles. 

For a network Q with conductances c on its edges, one can define the 
vi^eight w{T) of a spanning tree T to be the product of the conductances of 
its edges, w{T) = HeerCe- 

The most important theorem about spanning trees is that the Laplacian 
determinant is the weighted sum of spanning trees (this is also called the 
partition sum). 

Theorem 4.2 (Matrix- Tree Theorem [22]). Let deto A be the product of the 
nonzero eigenvalues (with multiplicity) of A for a connected network Q . Then 



-^detoA = 5^u;(T). 



Here the sum is over spanning trees of Q. 

This theorem is usually attributed to Kirchhoff [22]. An equivalent for- 
mulation is that the partition sum Z := X1t^(-^) ^^ equal to the reduced 
determinant of A, which is the determinant of A after removing any row 
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and column; since A is symmetric and constant functions are in the kernel 
any choice of row and column gives the same reduced determinant. 

Proof of Theorem \4.S^ sketch. Several proofs of the Matrix-tree Theorem are 
known; the classical one we present uses the Cauchy-Binet Theorem which 
is a formula for computing the determinant of a product of two nonsquare 
matrices: it states that if X is an n x m matrix and y is m x n with m > n 
then 

det XY = J2 det(X^) det(yB) (11) 

B 

where B runs over all choices of n column indices for X, X^ is the n x n 
minor of X defined by columns B and Yb is the nxn minor of Y defined by 
rows B. 

For the Matrix- Tree Theorem, let Aj j be the submatrix of A obtained 
by removing row 1 and column 1. Then Aj j = d^Cdi where di is obtained 
from d by removing column 1. 

One now uses X = d^ and Y = Cdi in fITT]) . Each B in flTTj) is a set of 
|y| — 1 edges; we claim that the nonzero summands in fITT]) are precisely the 
sets of edges which form spanning trees. If B forms a tree, its edges can be 
directed towards the removed vertex. Then in the definition of det dig as an 
expansion over the symmetric group, 

det dl^B = 2^ { — ^yd2a{2) ■ ■ ■ dna(n)-, 

the only nonzero term is the one whose permutation a matches each vertex 
to the edge adjacent to it and closer to the root. Thus det di^ = ±1. If edges 
in B do not form a tree, then there must be more than one component; the 
function on vertices which is 1 on a component not containing the removed 
vertex and zero elsewhere is in the kernel of dig, so det di^s = 0. 

Finally, the determinant of the diagonal submatrix det C^ for a tree B is 
the product of the edge weights, and det rfj ^ = det diB = ±1. So 

detidl^C^di^n) = w{T). 

This completes the proof. D 

The spanning tree measure is the probability measure on spanning 
trees of a graph in which the probability of a tree is proportional to its 
weight. In the special case that the conductances are all 1 this is the uniform 
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Figure 2: A grove of type 12|3 (nodes are in white; internal vertices in black). 



measure on spanning trees (UST measure); each tree has probability - 
where k is the determinant of the reduced Laplacian of the conductance- 1 
network; k is sometimes called the complexity of the network. 

Cay ley's formula [5], first proved by Borchhardt [3] gives the number 
spanning trees of the complete graph Kn to be n"~^. This follows easily from 
the above theorem since 



A 



(n-l -1 ... 
-1 n-l -1 



-1 \ 



-1 



\ -1 ... -1 n-l/ 
whose eigenvalues are n with multiplicity n — 1 and with multiplicity 1. 



4.3.2 Groves from minors of L 

On a network with nodes A^, a grove is a subset of the edges, each of whose 
components is a tree containing one or more nodes. In other words it is a 
"spanning forest" in which each component is connected to the boundary 
(and isolated vertices are considered components, so there are no isolated 
internal vertices). See Figure O We allow a component to contain more than 
one node; for example a spanning tree is a special case of a grove. 

Groves can be put into equivalence classes depending on how their com- 
ponents partition the nodes. Given a partition vr of the nodes, let St^ be the 
set of groves in which each component contains exactly all the nodes of a 
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part of vr. For example when there are three nodes the partitions are: 1|2|3 
where each node is in a separate component, 1|23 where nodes 2 and 3 are 
in the same component different from 1, and hkewise 2|13,3|12, and finally 
123 where all nodes are in the same component. More generally groves in 
the set 5*1121. ..|n; that is, in which each component contains exactly one node, 
are called uncrossings. The groves in Si2...n are spanning trees. 

Note that for circular planar graphs not all partitions vr are possible; for 
a partition to be obtained by a grove it must be planar. For example 13|24 
is not the partition of any grove on a circular planar graph. 

The probability measure on spanning trees extends to a probability mea- 
sure on groves: the weight of a grove is the product of its edge weights, and 
the probability of a grove is -| times its weight, where Z = '^tw{T) is the 
weighted sum of all groves. 

For each partition vr we denote by Z-j^ the weighted sum of groves of type 
TT. Then Z = YItt ^-^^ ^"^^ ^^e probability that a random grove has type tt is 
^. We define Zunc = Z\\2\,,\n to be the partition sum for uncrossings, and 
Ztrcc = ZYi...n to be the weight sum of spanning trees. 

All minors of the response matrix have combinatorial interpretations in 
terms of groves, for example when i ^ j we prove in Theorem 14.31 below that 

J _ ^7r(»i|rest different) _ Pr(7r(ij |rest different)) 
'^ Znnc Pr(unc) 

where 7r(ij|rest different) is the partition in which all nodes are in separate 
components except for i and j which are in the same component. More 
generally, we have the following combinatorial interpretation of all minors of 
L (a closely related result can be found in |i8j): 

Theorem 4.3 ([20]). Let Q,R,S,T C M be a partition of N into four sets 
with \R\ = \S\ (and some of which may be empty). Then det -^'1^^ is the ratio 
of two terms: the denominator is Zunc' the numerator is a signed weighted 
sum of groves of Qt, the graph Q in which all nodes in T are considered 
internal vertices, the nodes in Q are in singleton parts, and in which nodes 
in R are paired with nodes in S, with the sign being the sign of the pairing 
permutation: 

7 rSp(i) I |*p{fc)| I I 1 
detL|- = (-l)m Y: (-ir ^" zj ■ ^''^ 

permutations p 
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As a special case of this theorem, we have the above statement about Lij 
for i ^ j (here L^j is the i, j-entry of L which in the notation of the theorem 
would be L\'^\)- Another special case is the reduced determinant of L: for 
any i and j, removing row i and column j gives 

detLU^ 



(to see this, take R = {i}, S = {j} and T = Af \ {i,j}). 

In the example (jH]) above where Q is the Y network of Figure [H L12 



ClC2 



, , -; C1C2 is the weight of the unique grove of type 12 1 3 and ci + C2 + C3 
is the sum of the weights of the three groves of type 1|2|3. 



Proof of Theorem\4.1\ necessity. Since nonplanar pairings do not occur, if R 



and S are non-interlaced subsets of the nodes of the same size, then det L;| 
is a nonnegative because there is a single term in the sum in flT2|) and the 
permutation is the identity. This gives a combinatorial interpretation of the 
inequalities in Theorem 14. 11 they are positive because they count something. 

D 

Proof of Theorem \4.3\ Write A in block form as 

where A is indexed by the nodes and C by the internal vertices. Then we 
can factor A as 

A _ M ^\ - (A-BC-^B* BC-^\ (I 0\ _ (-L BC-^\ (I 
^~\B' C)~\ I )\B' C)~\Q I )\B' C 

We wish to compute the determinant of the minor A^^^. For notational 
simplicity reorder the nodes so that nodes in X come last and nodes in Y 
come first. Then 



A 

and 

^yu7 _ / ~Lx * 




^XUI 








The Cauchy-Binet formula fill I) then gives 



In particular 



det A^^^^ = (-1)1^1 det L^ det C. 



det Al^^^i = (-l)l«l+l^l det Ll\^l det C 



where / is the set of internal nodes. Since Zunc = det C, it suffices to show 
that det A^^\^^j is the desired numerator. 

Let J' = / U T be the "new" set of internal nodes. We need to evaluate 
detA'^^^,. The remainder of the proof is now an extension of the proof 
of Theorem lO Write A = d*Cd. Then A|^^', = d*s^jj,CdRui' where dx 
represents the restriction of d to the subspace indexed by X. By the Cauchy- 
Binet theorem, 

det d*sijj,CdRijr = 2_^ det((i|^jj/)* det Cy det d^ijj,, (13) 

Y 

where the sum is over collections of edges Y of cardinality ISU/'I = |-RU/'|. 
The terms in the sum f lT3|) are collections of edges Y in the graph Q. The 
number of components of the subgraph with edges Y is at least 

|vertices| - |edges| = \RU S U Q U I'\ - \Y\ = \R\ + \Q\. 

We claim that for each nonzero term in the sum f ll3p . each component is a 
tree containing exactly one element oi S U Q. If some component did not 
contain a point of 5* U Q, then the function which is 1 on this component 
and zero elsewhere would be in the kernel of d^^jj, whose determinant would 
then be zero. Since the number of components is \R\ + \Q\ = \S\ + \Q\ 
each component must contain exactly one point of S U Q. Similarly each 
component also contains exactly one element oi RU Q (since detd'^^jj ^ 0). 
Thus each component is a tree containing either a unique q & Q (and no 
point of -R U S"), or a tree containing both a unique r E R and a unique 
s E S. The total number of components is \R\ + \Q\ so each element of Q, R 
and S* is in a component. In particular there is a permutation p and for each 
i a component joining r^, the ith element of R, to Sp(j), the p(z)th element 
of S. The weight of a component is the product of its edge weights. The 
sign of the term corresponding to the subset Y can be determined as follows. 
Each component connected to Q occurs in both det d^yjp and det d^yjj, and 
so contributes sign +1. If we postmultiply dsur by the permutation matrix 
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P which permutes the rows according to p then the sign of the corresponding 
terms in Pdsur and djujjr are the same since they corresponding to the same 
choice of rows. Since detP = (—1)^ this completes the proof. D 

4.3.3 Other groves 

Theorem 14.31 shows that, when divided by the factor Pr(unc) (the probabil- 
ity of the uncrossing) certain grove probabilities are given by noninterlaced 
minors of L: those groves which correspond to noninterlaced pairings, for 
example. 

Remarkably, one can give a similar expression for all grove types: see 
Theorem 14.41 below. 

For a partition r of {1, . . . , n} (planar or not) we define 

^-=E n ^M' (14) 

F {i,j}GF 

where the sum is over those spanning forests F of the complete graph on n 
vertices 1, . . . ,n for which the trees of F span the parts of r. For example 
we have 

-^(123|4) = L12L13 + L12L23 + L13L23 

and 

-^(13|24) = L13L24. 

Theorem 4.4 ([19]). For any planar partition a we have 

Pr((T) 



Pr(unc) 



J^^a.r^r 



where V is the integer matrix defined below. 

The rows of the matrix V are indexed by planar partitions, and the 
columns are indexed by all partitions. In the case of n = 4 nodes, the 
matrix V is 
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CO 



^^c^^co Oi^ CO ^ ^ ^ CO 

COCMCM^^^COCMCOCOCMCM 



CM CM CO 



CM CO ^ 



CO ^:-!- 
CM CO 



1|2|3|4 10000000000000 

12|3|4 01000000000000 

13|2|4 00100000000000 

14|2|3 00010000000000 

23|1|4 00001000000000 

24|1|3 00000100000000 

34|1|2 00000010000000 

12|34 00000001000000-1 

14|23 00000000100000-1 

1|234 00000000010000 

2|134 00000000001000 

3|124 00000000000100 

4|123 00000000000010 

1234 00000000000001 

For example, the row for 1|234 tells us 
Pr(l|234) 



1 
1 
1 
1 
0. 



Pr(l|2|3|4) 
and the row for 12134 tells us 



-^(l|234) + -^(13|24) — -^2,3-^3,4 + -^2,3-^2,4 + -^2,4-^3,4 + -^1,3-^2,4 



Pr(12|34) 



L 



(12|34) 



L 



(13|24) 



-^1,2-^3,4 — Li^L 



2,4- 



Pr(l|2|3|4) 

We call this matrix V the projection matrix from partitions to pla- 
nar partitions, since it can be interpreted as a map from the vector space 
whose basis vectors are indexed by all partitions to the vector space whose 
basis vectors are indexed by planar partitions, and the map is the identity 
on planar partitions. For example, the column for 13|24 tells us 

13|24 projects to - 12|34 - 14|23 + 1|234 + 2|134 + 3|124 + 4|123. 

The projection matrix may be computed using some simple combinatorial 
transformations of partitions. Given a partition r, the rth column of V may 
be computed by repeated application of the following transformation rule, 
until the resulting formal linear combination of partitions only involves planar 
partitions. The rule is a generalization of the transformation 

13124 -> -12134 - 14123 + 11234 + 21134 + 31124 + 41123 
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which can be represented as a sort of skein diagram: 




+ 





+ \ + / + 





We generahze this rule to partitions r containing additional items and parts 
as follows. If partition r is nonplanar, then there will exist items a < b < 
c < d such that a and c belong to one part, and b and d belong to another 
part. Arbitrarily subdivide the part containing a and c into two sets A and 
C such that a & A and c & C, and similarly subdivide the part containing b 
and d into B 3 b and D 3 d. Let the remaining parts of partition r (if any) 
be denoted by "rest." Then the transformation rule is 

AC\BD\ rest -^ A\BCD\ rest + B\ACD\ rest + C\ABD\ rest + D\ABC\ rest 

- AB\CD\ rest - AD\BC\ rest . (Rule 1) 

Theorem 4.5 ([I9]). Any partition r may be transformed into a formal 
linear combination of planar partitions by repeated application of \Rule l\ and 
the resulting linear combination does not depend on the choices made when 
applying \Rule 1[ so that we may write 

planar partitions a 

4.4 Minimality and electrical equivalence 

4.4.1 Medial graph 

If ^ is a circular planar network or network on a surface, the medial graph 
of Q is the graph with a vertex t>(e) for every edge e of Q, and an edge 
connecting v{e) and v{e') if e and e' share a vertex and are adjacent edges 
(in cyclic order) out of this vertex. The medial graph is regular of degree 4. 
We usually cut the edges of the medial graph which separate the nodes from 
oo, so that the medial graph has two half-edges (called "stubs") adjacent to 
each node, see Figure [31 

A strand of the medial graph M is a maximal path in the medial graph 
which continues "straight" at each vertex of M, that is, turns neither left nor 

22 



@K^-— © 



© '^-;C_^__ i J><i_ i J::-<' © 



Figure 3: The medial graph (dotted) for a 4-node network (sohd). 

right. There is a strand starting at each stub and ending at another stub, 
and possibly other strands which form closed loops. 

If Q is drawn in a disk, with the nodes on the bounding circle, the strands 
of a medial graph are the boundaries of a cell decomposition of the disk. The 
cells come in two types, corresponding to vertices and faces of Q (between 
two adjacent nodes there is a cell corresponding to an exterior face of Q). 

A nice property of medial strands is that, for minimal graphs (see the 
next section), they are boundaries of the support of equilibrium potentials. 
Let us explain. 

Lemma 4.6. JE/ Suppose Q is minimal. Let a, h he stubs connected by a 
medial strand. This medial strand divides the disk into two parts U+,U^. 
Then there is an equilibrium potential f (a function harmonic on I and with 
Dirichlet boundary conditions on M), not identically zero, with harmonic 
conjugate f* with the property that both f and f* are identically zero on 
(cells in) f/_. 

For an example see Figure IH Furthermore we have 

Lemma 4.7. JE/ If G is minimal and a, b are two stubs, let C+ and C_ be 
the two arcs of the boundary of the disk separated by a and b. Then a and 
b are connected by a medial strand in Q if and only if there are equilibrium 
potentials /+ and /_ and conjugates respectively g+,g- such that f+,g+ are 
zero on nodes and exterior faces of C^, and f-,g- are zero on nodes and 
exterior faces of C-. 
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Figure 4: A network (solid) with a strand (dashed) and an equihbrium po- 
tential and dual which are zero on one side. In this example all conductances 
are 1. 

4.4.2 Minimality and electrical transformations 

A circular planar network is said to be minimal if the medial strands have 
the properties 

1. There are no closed loops. 

2. A strand does not intersect itself. 

3. No two strands intersect more than once. 

The reason we define minimal graphs is that every circular planar network 
is equivalent, in the sense of having the same response matrix, to a minimal 
circular planar network. Moreover minimal graphs are precisely those in 
which the reconstruction has a unique solution [6]. 

In [7] it was shown that any circular planar network can be converted 
into a minimal circular planar network by a sequence of local transformations 
called electrical transformations. An electrical transformation is a local 
rearrangement of the graph and conductances of one of the types shown in 
Figure [51 They consist of 

1. Removing a "dead branch" (remove a non-node vertex of degree 1 and 
the edge connecting it to the rest of the graph) or a self-loop. 
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Figure 5: The electrical transformations. 



2. Replacing two conductors ci,C2 in series with a single conductor of 



conductance 



C1C2 
C1+C2 



(as long as the common vertex is not a node). 



3. Replacing parallel edges of conductances Ci,C2 with a single edge of 
conductance ci + C2. 

4. the "Y-Delta" move, also called star-triangle move: an internal vertex 
of degree three with edges of conductance Ci,C2,C3 is replaced by a 
triangle with edges of conductance — r^, — , — r^, — , — r^, — as in the 

O O C1+C2+C3 ' C1+C2+C3 ' C1+C2+C3 

diagram. 

The reverse of a Y-Delta move is a Delta-Y move and is also allowed. The 
inverse transformation takes conductances ci,C2,C3 on the "Delta" to con- 
ductances 

C1C2 + C1C3 + C2C3 C1C2 + C1C3 + C2C3 C1C2 + C1C3 + C2C3 



Cl 



C2 



C3 



on the Y. For simplicity when we say "Y-Delta move" we mean either a 
Y-Delta or its reverse. 

We say that circular planar networks Q,Q' on the same number of nodes 
are topologically equivalent if one can be converted to the other as graphs 
(that is, ignoring conductances), using electrical transformations. We define 
circular planar networks Q, Q' to be electrically equivalent if they have 
the same response matrix. 

4.4.3 Topological equivalence 

It is not hard to show that any circular planar network is topologically equiv- 
alent to a minimal circular planar network, that is, can be converted to a 
minimal network by electrical transformations. The idea is to simply chart 
how the moves act on medial strands: Each of the first three types of moves 



25 



decreases the number of crossings of the medial strands (that is, the num- 
ber of edges of the network). We isotope the strands, leaving their endpoints 
fixed, until each strand has no self-intersections and crosses each other strand 
at most once. During this isotopy, each time a strand crosses an intersection 
of other strands we perform a Y-Delta move on the graph; other singular- 
ities encountered will decrease the number of crossings and so are finite in 
number. 

This argument shows that for the purposes of determining topological 
equivalence one can assume Q and Q' are minimal. What is a bit harder to 
show is that the "minimization" process leads to a unique minimal graph, up 
to Y-Delta moves. That is, if Q' and Q" are minimal graphs obtained from 
the same graph Q then Q' can be converted to Q" using only Y-Delta moves 
(and Delta-Y moves). See Theorem 14.81 below. 

For a given minimal circular planar graph Q with n nodes, the medial 
graph has 2n stubs, one starting to the left and right of each node. The 
medial paths connect these stubs in pairs. Let vr denote the pairing of the 
stubs by the medial paths; if we label the stubs from 1 to 2n in cclw order 
then we can think of vr as a fixed-point free involution of {1, 2, ... , 2n}. 

Theorem 4.8 ([6J). Two minimal graphs Q,Q' are topologically equivalent if 
and only if they have the same medial strand stub involution. 

Proof. If two minimal graphs have the same medial strand involution then 
one can isotope the strands without forming any new crossings so the two 
graphs are isomorphic: for example put the nodes and stubs in general po- 
sition on a circle respecting their circular order; now isotope the strands to 
straight chords across the disk. 

Conversely, by Lemma 14.71 the response matrix of a minimal graph Q 
determines its stub involution, since the supports of equilibrium potentials 
and their harmonic conjugates are determined by L. D 

4.4.4 Well-connected and non-well-connected networks 

Recall that a network is well-connected if all non-interlaced minors of L are 
strictly positive. Equivalently, for any non-interlaced subsets of nodes A, B, 
there is a grove connecting points of A to those of B in pairs. (If |^| = \B\ = k 
one can just look for a set of k pairwise disjoint paths connecting points of 
A to those of B] such a set is easily completed to a grove by adding extra 
edges.) 
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Figure 6: The standard networks Tk for k 
connected. 



6. These are well- 







Figure 7: Some (almost) circularly symmetric well-connected graphs. 

What stub involution corresponds to a well-connected network? It is not 
hard to construct well-connected graphs. One nice family is illustrated in 
Figure [61 We leave it to the reader to show that these graphs are well- 
connected, and have stub involution 7r„eii which pairs i with n + i for each i. 
That is, each stub is paired with the diametrically opposite stub. 

Corollary 4.9. For a minimal graph Q on n nodes, the following are equiv- 
alent: 

1. Q is well-connected. 

2. Its stub involution is 'Kweii- 

3. Q is topologically equivalent to r„. 

Another nice family of well-connected graphs which are in addition cir- 
cularly symmetric (for n odd) and nearly circularly symmetric (for n even) 
is given in Figure [71 See Figure [HI for a third family. 

What about non-well connected graphs? Is there any structure to the set 
of topological equivalence classes of n-node networks? See below. 

4.4.5 Electrical equivalence 

We discussed topological equivalence above. Regarding electrical equivalence, 
we have 
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Figure 8: A different family of well-connected graphs for even n. 

Theorem 4.10. Electrical transformations do not change the response ma- 
trix. 

Proof. It suffices to show that electrical transformations preserve current 
flow. More precisely, the current flow on the "before" network will be equal 
to the current flow on the "after" network if defined appropriately on any 
added edges. 

This is clear for transformations of type 1 (removing a dead branch or 
self-loop) since these edges have no current fiow. For a type 3 move (parallel 
reduction) the current fiow on ci is Cidf{e) where / is the corresponding 
potential on the vertices. Similarly the current fiow across C2 is C2df{e). The 
current fiow on the combined edge is (ci + C2)df{e) which is the sum of the 
current fiows on the individual edges. 

For a type 2 (series) transformation, the relevant equation is Cidf{ei) = 
C2(i/(e2) since no current is lost at the center vertex. This implies that 



cidf{ei) 



_ClC2_ 
Cl +C2 



(c?/(ei) + rf/(e2)). 



and df{ei) + df{e2) is the potential drop across the combined edges. 

A similar computation works for the Y-Delta move. When converting 
from a Y to a Delta, the current along the edges of the Y directed towards 
the central vertex is Cidfi, C2df2, c^df^ (where the dfi are the potential drops) 
which must add to zero since no current is lost at the central vertex: 

Cidfi + C2df2 + csdfs = 0. 



28 



A short computation then shows that 

cidfi = Ci2idfi-df2) + Ci3idfi-dfs), 
C2df2 = Cuidh-df2) + C,sidh~dfs), 
czdh = Cuidf, - df2) + C,,idf, - df,) 

for arbitrary dfi, df2, df2 when and only when Cij = ^ 1'^^^. ■ D 

As a consequence of this theorem, for determining the space of response 
matrices of planar networks we need only consider minimal networks. 

4.5 Reconstruction and the space of networks 

Given a matrix which satisfies the conditions and inequalities of Theorem 
(14. ip . is it the response matrix of a circular planar network? How can one 
determine the conductances from L? This is the Reconstruction Prob- 
lem or EIT problem. It is closely related to the Electrical Equivalence 
Problem: what sets of minimal networks are electrically equivalent? 

4.5.1 Algorithm 

Both [8] and [6] gave an iterative algorithm for reconstruction on any minimal 
graph. The idea is as follows. Any network is equivalent to one which has 
either an edge connecting adjacent nodes or a node of degree 1 (take a strand 
ab with the property that there are no strands completely contained on one 
side f/+ of it; using Y-Delta moves, move all crossings in U+ of other strands 
to the other side f/_. The node adjacent to a and inside U+ will either have 
degree 1 or be connected to its adjacent node on the other side of a). 

One can now compute the conductance on this edge e from L: find the 
potential with values zero on side U^ of strand ab. Then the conductance Cg 
is easily obtained from the current. 

Removing the edge e (if it connected nodes) or contracting the edge e 
(if it is adjacent to a degree- 1 node) results in a new network which is still 
minimal, and the new response matrix is a simple rational function of the 
old response matrix. 

This algorithm shows that the conductances are rational functions of the 
Lij. However it is not so easy to use this to give an explicit formula for the 
conductances as a function of the response matrix. In [20J an explicit formula 
is given for the reconstruction problem for the standard graphs F^ of Figure 
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El which expresses Cg as a ratio of polynomials in the Lij which are Pfaffians 
of skew-symmetric matrices formed from L. This can be considered in some 
sense the best possible result for r„ since the polynomials of the Lij are 
shown there to be irreducible in general. However this does not preclude the 
possibility that for other well-connected networks there is a simpler formula: 
see 



4.5.2 Partial order on networks 

By allowing some conductances to go to zero and others to tend to oo a 
network degenerates (topologically, we remove edges of conductance zero 
and contract edges of conductance oo). 

There is a partial order on networks under this degeneration. We define 
Q < Q' iia. network equivalent to Q can be obtained from a network equivalent 
to Q' by contraction and deletion of edges. 

One can show that every minimal network is a minor (in the sense of 
contraction and/or deletion of edges) of a well-connected network. So the 
well-connected network is the unique maximal element in this partial order. 
If we allow disconnected networks, the networks which are minimal are those 
whose stub involution is a non-crossing matching; see Figure [TOl 

Lemma 4.11. If Q,Q' are minimal circular planar networks and Q < Q' 
then there is a sequence Qi < Q2 ^ ■ ■ ■ "^ Qk = Q' of minimal networks where 
Qi is equivalent to Q and Qi is obtained from a network equivalent to Qi^i by 
deletion or contraction of a single edge. 

Proof. First suppose that Q is obtained from Q' by deleting a single edge e 
(the case of contraction is equivalent by duality). Here is how one minimizes 
the resulting network: let a, 6, c, d in cclw order be the stubs of the strands 
ac and hd crossing at edge e. Let t/i, t/2, ^3, U^ be the four arcs of the circle 
separated by a, 6, c, (i, so that Ui is the arc from a to 6, etc. Using Y-Delta 
moves we can first comb strands so that strands connecting stubs from Ui 
to U^ cross "left" of e, that is, do not enter the region delimited by the 
path 6ec, and strands connecting stubs from U2 and t/4 cross "below" e, that 
is, do not enter the region delimited by the path ced. See the Figure [91 
Now we resolve the crossing at e, creating a strand connecting a and h and 
another connecting c and d. The strand ah forms a bigon with any strand 
from f/2 to f/4. To minimize the resulting network we resolve all crossings 
of these strands with the semi-strand from 6 to e as indicated in Figure [HI 
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The resulting network is now minimized (and any minimization results in a 
network equivalent to this one). If we perform this sequence of resolutions 





Figure 9: When a conductance e goes to zero, comb strands left and down 
(left figure); then resolve strand crossings between b and e (right figure). 

in order of the crossings along the semistrand be from the boundary towards 
the center, (resolving the crossing at e last) the intermediate networks are 
all minimal. This completes the proof in the case when Q is obtained from 
Q' by resolution of a single crossing. 

In the general case, we contract and delete multiple edges of Q', then 
minimize to get Q. Take a path in the space of conductances which sends the 
desired conductances to zero or infinity, one after another, in any order. Let 
L[t) be the corresponding path in the space of response matrices. When we 
minimize the resolution of a single crossing as above, the resulting L matrix 
does not change (in the sense that no additional minors become zero), since 
minimization leaves L invariant. Thus the sequence of resolutions is replaced 
with a sequence of "minimal" resolutions, that is, resolutions which preserve 
minimality. D 

This partial order on networks induces a partial order on fixed-point free 
involutions of {!,... ,2n}. This partial order is graded by the number of 
crossings in the involution, that is, the number of pairs i < j < k < i with 
i paired with k and j paired with i. Moving one level down in the partial 
order corresponds to resolving exactly one crossing which is adjacent to the 
boundary (in the sense that there is a path in the disk from the crossing to 
the boundary circle which does not meet any strand). Algebraically, moving 
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up one level corresponds to the operation of conjugating the involution by a 
transposition in such a way as to increase the number of crossings by 1. 
See Figure [TO]for the Hasse diagram of the partial order in the case n = 3. 




Figure 10: The partial order on circular planar networks with three nodes. 
Only the medial graphs are shown. 

This partial order describes the cell structure of the space of networks. 

4.5.3 Minimal sets of inequalities 

What is the structure of the set of response matrices of all circular planar 
networks? The set Q of response matrices of circular planar networks on n 
nodes is a closed subset of M'^("~i)/2 if we use the coordinates {-^ij}rt>i>i>i- 
It is a semi- algebraic set, that is, defined by a set of algebraic inequalities 
det L^ > for the non-interlaced subsets A, B of nodes. Its interior Vt^ is 
defined by using strict inequalities, and corresponds to the well-connected 
networks. 
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The number of inequalities defining Q is exponential in n. However there 
is a smaller set of n{n — l)/2 inequalities which defines Q'^ . See below. 

If we consider a topological equivalence class (defined by a stub invo- 
lution tt) of graphs, the response matrices of networks supported on these 
graphs form a subset Q-j^ G Q, which is a subset of the boundary of fl if 
71 ^ 7r„cii- fi has the structure of a cell-complex, in which the Q-,^ are the 
cells, see Lam and Pylyavskyy [23j . Is each cell of dimension k defined by k 
minor inequalities (and n{n — l)/2 — k equalities)? The answer is probably 
yes, but at the moment our understanding is limited. A similar situation 
where the cell structure is explicitly worked out is Postnikov [29] who deals 
with totally positive/totally nonnegative matrices and the totally nonnega- 
tive Grassmannian. 

Let us discuss here only the minimal sets of inequalities defining the 

set r2+. Somewhat remarkably, there are many different sets j J minor 

determinants of L whose positivity implies the positivity of all non-interlaced 
minors. The situation resembles that of a cluster algebra but so far we have 
been unable to find the relevant cluster structure. 

One easy-to-remember set of I „ ) minors are the central minors, de- 
fined as follows. Suppose first n is odd. Define Mj^i for 1 < i < {n — l)/2 to 
be the minor L^ with A = {1,2, . . . ,i} and B = {{n - l)/2 + 1, (n - l)/2 + 
2, . . . , (n — l)/2 -|- i}. Let Mj ,,■ be the minor obtained from Mj i by rotating 
the indices cyclically by j — 1, so that A = {j,j + 1, . . . ,j + i} and likewise 
for B. 

These minors are called central minors since, if we arrange the n points 
evenly spaced on a circle, and connect each element of A to the corresponding 
point of B opposite it, we get a set of i parallel chords which are "central" 
in the sense that they are as close to a fixed diagonal of the circle as possible 
while remaining disjoint from each other. See Figure [TT] for an example. In 
the case n is even we modify the above; define Mij for i even as before, 
but for i odd we have, for each diagonal, two choices of positioning of the 
chords to make them "most central" ; choose one of these arbitrarily. These 
"off-center" minors are part of our set of central minors. See Figure [T^] for a 
natural choice. 

Theorem 4.12 ([2I])- If the central minors {Mij} are positive then all non- 
interlaced minors are positive. 
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Figure 11: The 21 central minors for a well-connected network on 7 nodes. 
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Figure 12: The 15 central minors for a well-connected network on 6 nodes. 
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The following proof is lifted directly from 



Proof. This follows using two identities for minors. 

Let [/ be a matrix, and let a, b, c index some of its columns, and z, d index 
some of its rows. Then it is elementary that 



'^ z,a '^ z,b '^ z,c 
^z,a '^ z,b '^ z,c 
Ud,a Udfi Ud^c 



U.^a 



Uz,h Uz^c 

Ud,b Ud,c 



U: 



z,b 



'-^ z,a '^ z,c 
Ud,a Udr 



+ u. 



Uz,a Uz^b 
Ud,a Udh 



Suppose U is invertible, and let M denote its inverse. Let ri, . . . ,rfc denote 
the row indices of M (and column indices of [/), and ci, . . . , c^ denote the 
column indices of M (and row indices of U). Dividing through by (det?7)^ 
and using Jacobi's identity relating minors of a matrix to minors of its inverse, 
we obtain 

= det(M,,);5: det(M,,);5;,-det(M,,)i5^ det(M,,)g:;^+det(M,,);5: det(M,,)55:^ 

Since the column z is always excluded, we may as well suppose that M 
is a fc X (A; — 1) matrix. Then 



^^^6 



\ij^a,c 



\iy^a 



iy^b,c 



\ij^c 



\i^a,b 



det(M,,)r" det(M,,)5r = det(M,,)f'^ det(M,,)5r+det(M,,)f ^ det(M,,);y/- 
For example, taking M = i^i'2'3'4' and a, b,c,d = 9, 8, 5, 4, we have 



7-9,7,6,5 

1,2,3,4 




7-8,7,6,5 
1,2,3,4 




7-9,7,6 
1,2,3 


+ 


7-9,8,7,6 
1,2,3,4 




7-7,6,5 
1,2,3 










r8,7, 
^1,2, 


5 
3 









We can denote this pictorially as 



m 



We call transformations of this type the "jaw move". 

For a given non-interleaved determinant interspersed with at least one 
isolated node, we can take b to be one of the interspersed isolated nodes, and 
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a and c to be the first and last of tlie nodes on tlie same side as 6, and d to be 
eitlier first or last node on the other side as h. With this choice of a, 6, c, d, the 
jaw move expresses the original determinant as a positive rational function of 
"simpler" non-interleaved determinants, where a determinant is simpler if it 
has fewer strands, or else the same number of strands but fewer interspersed 
isolated nodes. 

By repeated application of the jaw rule, any non-interleaved determinant 
can be expressed as a positive rational function of non-interleaved contiguous 
determinants. 

The other identity that we need is based on Dodgson condensation (and 
is also called the Desnanot-Jacobi identity). It can be derived in a similar 
way as the jaw move, with the elementary starting identity 



UrM, 



ca^dM 



Uc,a Uc,b 
Uda Udb 



+ UchUd,a- 



Dividing through by (det f/)^ and using Jacobi's identity, we obtain 
det(M,,,)55:det(M,,)g^ = det(M,,,)} det(M,,);g+det(M,,)gtdet(M,,,)g:. 
For example, taking M = L^^z ^"^^ '^) b,c,d = 8, 6, 3, 1, we obtain 



r8,7 




7-8,7,6 
-^1,2,3 




LI 


+ 


r7,6 

-^^1,2 




r8,7 
-^2,3 










r7,6 
-'^2,3 









which we can denote pictorially as 




We call transformations of this type the "condensation move" . 

If both sides of the determinant are contiguous, but the crossings are 
off-center, then we can pick a and d to be the paired nodes that are most 
off-center, and b and c to be the nodes next to the crossing but not in the 
crossing (and closer to the center). Then after a condensation move, the 
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original determinant is expressed as a positive rational function of contiguous 
crossings that are strictly more central. By repeated application of such 
condensation moves, any contiguous determination may be expressed as a 
positive rational function of central contiguous crossings. 

The remaining case to check is for n even, with the (minimally) off-center 
contiguous crossings come in pairs, of which only one is included in the base 
cluster. Using a condensation move we can express a minimally off-center 
contiguous crossing in terms of its opposite minimally off-center contiguous 
crossing and central contiguous crossings as shown below 




n 



4.6 The Jacob ian 

A remarkable property of the birational map from conductances to the (ap- 
propriate set of) L matrix minors for minimal networks is that the Jacobian 
is ±1: 



Theorem 4.13 ([21]). 



O log Cj 



In other words the volume form Y\^ — on the space of conductances is 
mapped to ±1 times the volume form Y\_i -jj^ on the space of L-matrix minors. 
Note that the matrix entries have a probabilistic interpretation: each Mj can 
be written Mj = Pr(7r)/Pr(unc) for some noninterlaced pairing vr. Thus 



dlog Mi 
d log Cj 



Cj 



dZ{ 



vr 



dZ,, 



Z(n) dc, 



dci 



is the difference in the probability that edge Cj is in a random vr-pairing, 
minus the probability that Cj is in a random uncrossing. 
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The proof of this theorem rehes of techniques developed in [22] and 
and is too long to give here. Let us just give an example, for a well-connected 
network on 3 nodes, the Y of Figure [H The central minors in terms of the 
conductances are La = — p-^, — . The above Jacobian matrix (with columns 

'J C1+C2+C3 \ 

in the order L12, L13, L23) is 



C1+C2+C3 C1+C2+C3 C1+C2+C3 

1 1 1 1 



CI+C2+C3 CI+C2+C3 Cl+C^+C3 

1 1 1 



1 



C1+C2+C3 C1+C2+C3 C1+C2+C3/ 

whose determinant is 1. 

5 Networks on surfaces with nontrivial topol- 
ogy 

Some of the material in the preceding section on circular planar networks has 
been extended to the case of annular networks and more generally networks 
on surfaces. However the theory is not as complete at present as in the case 
of circular planar graphs. Moreover beyond the annulus and torus the theory 
is for the most part nonexistent. 

We can nonetheless ask the same questions as arose in the planar case: 

1. What is the natural notion of response matrix? 

2. For which networks can we reconstruct the conductances from the re- 
sponse matrix? 

3. What is the structure of the space of response matrices? 

4. What natural combinatorial objects are relevant? 

Our answers, briefly, are as follows. We elaborate on these in the following 
sections. 

1. We use (the Schur reduction of) the bundle Laplacian for a fiat C*- 
connection or a fiat SL2(C)-connection on the network. 

2. We can conjecturally reconstruct conductances for minimal networks 
on the annulus. For other surfaces reconstruction is not in general 
possible, unless we are given more data: we (conjecturally) need to 
provide "spectral" data as well. 
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3. We don't know what the structure of the space of response matrices is, 
even for the annulus. However we have a description in the case when 
there are no nodes, for the annulus and torus. 

4. The natural combinatorial objects are "cycle-rooted spanning forests": 
forests in which each component has a unique (topologically nontrivial) 
cycle. 

5.1 Vector bundle Laplacian 

In the case of a network on a surface it is natural to consider not the standard 
Laplacian but a more general Laplacian operator, the vector bundle Lapla- 
cian, which depends on a connection on the bundle (see definitions below). 
This is because the standard Laplacian does not contain very much informa- 
tion about the topology of the surface on which the network is embedded... in 
particular the reconstruction problem is not solvable in general (as an exam- 
ple consider the network on an annulus with two vertices, both nodes, and 
two edges joining them with different homotopy classes). From the bundle 
Laplacian we construct a response matrix which contains more information 
than the one constructed from the standard Laplacian, and conjecturally 
allows reconstruction for minimal networks. 

Let us first define vector bundles and connections on a network. 

5.1.1 Vector bundles and connections 

Given a fixed vector space W, a W^-bundle, or simply a vector bundle 
on a network Q is the choice of a vector space W^ isomorphic to W for 
every vertex v oi Q. A vector bundle can be identified with the vector space 
Wg := ®vW,u = W^l^l, called the total space of the bundle. A section of a 
vector bundle is an element of Wg. 

A connection $ on a W^-bundle is the choice for each edge e = vv' 
of Q of an isomorphism 0^,^,/ between the corresponding vector spaces 0^' : 
W^ — )■ W^i, with the property that 0^,^,/ = 0~,J,. This isomorphism is called the 
parallel transport of vectors in W^ to vectors in W^i. Two connections $, <!>' 
are said to be gauge equivalent if there is for each vertex an isomorphism 
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ipv : Wy — )■ W^ such that the diagram 



W^ ^^^^^ W^. 



tpv 



^v' 



commutes. In other words $' is just a base change of $. Note that the 
connection has nothing to do with the conductances. 

Given an oriented cycle 7 in ^ starting at v, the monodromy of the 
connection around 7 is the element of End(W^„) which is the product of the 
parallel transports around 7. Monodromies starting at different vertices on 
7 are conjugate, as are monodromies of gauge-equivalent connections. 

A line bundle is a l^-bundle where W = C, the 1-dimensional complex 
vector space. In this case given a connection if we choose a basis for each C 
then the parallel transport along an edge is just multiplication by an element 
of C* = C \ {0}. The monodromy of a cycle is in C* and does not depend on 
the starting vertex of the cycle (or gauge). 

5.1.2 Flat connections 

A connection on a network embedded on a surface E is said to be a flat 
connection if the monodromy around faces of Q is trivial. This implies that 
the monodromy is trivial around loops on Q which are nuU-homotopic as 
loops on S. 

It is not hard to see that flat connections, modulo gauge equivalence, are 
in bijection with homomorphisms of vri(S) into the group of automorphisms 
of W^. 

5.1.3 The Laplacian 

The Laplacian A on a W^-bundle with connection $ is the linear operator 
A:Wg^Wg defined by 

where the sum is over neighbors v' oi v. 
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Note that if the vector bundle is trivial, in the sense that 0^^,/ is the 
identity for all edges, this is our the notion of graph Laplacian from ([T]) (or 
more precisely, the direct sum of dim W copies of the Laplacian) . 

As in the case of the standard network Laplacian we often think of A as a 
matrix A = (A^^/)^,^^,/^^ whose entries are —c^v' times the parallel transport 
from f ' to f . In particular the conductance is acting as a scalar multiplication 
inW. 

Here is an example. Let Q = K3 the triangle with vertices {f 1, f 2, "ya}- Let 
$ be the line bundle connection with (pv.v = Zij G C*. Then in the natural 
basis, A has matrix 



(C12 + Ci3 —C12Z12 — Ci32;i3 
-Cl2^r2^ C12 + C23 -023^23 | • (15) 

-CnZ^^ -023^23^ Ci3 + C23, 

5.1.4 Edge bundle 

One can extend the definition of a vector bundle to the edges of ^. In this 
case there is a vector space We — W for each edge e as well as each vertex. 
One defines connection isomorphisms (p^e = (pev ^^^ ^ vertex v and edge e 
containing that vertex, in such a way that if e = vv' then 0„„/ = 06^' ° 4>ve, 
where 0„„' is the connection on the vertex bundle. 

The vertex/edge bundle can be identified with lyl^l+l^l = Wg ® We, 
where We is the direct sum of the edge vector spaces. 

A 1-form (or cochain) is a function on oriented edges which is antisym- 
metric under changing orientation. If we fix an orientation for each edge, 
a 1-form is a section of the edge bundle, that is, an element of W^^^. We 
denote by A^{Q, $) the space of 1-forms and A^{Q, $) the space of 0-forms, 
that is, sections of the vertex bundle. 

We define a map ci : A^iG, $) ^ A\g, $) by df{e) = <PyJ{y) - 0,e/(x) 
where e = xy is an oriented edge from vertex x to vertex y. We also define 
an operator d* : A^ — )■ A*^ as follows: 



d*uj{v) = J2 



,,uj{e 



where the sum is over edges containing v and oriented towards v. Despite the 
notation, this operator d* is not a standard adjoint of d unless (pev and (p^e 
are adjoints themselves, that is, if parallel transports are unitary operators 
(see below). 

41 



The Laplacian A on A° can then be defined as the operator A = d*Cd as 
before: 

d*Cdf(v) = ^ Ce(t)evdf{e) 

e=v'v 
= ^ Ce(l)evi(f)vef{v) " (f^v'efiv')) 



e=v'v 
v' 



= A/(t;). 

We can see from the example (ITB]1 above on K^ that A is not necessarily 
self-adjoint. However if 0^^/ is unitary: 0~J; = 0*,^ then d* will be the adjoint 
of d for the standard Hermitian inner products on W'^' and W'^', and so in 
this case A is a Hermitian, positive semidefinite operator. In particular on a 
line bundle if \(f)vv'\ = 1 for all edges e = vv' then A is Hermitian and positive 
semidefinite. 

5.2 Minimal networks 

Here we show that any network on a surface can be "minimized", that is, 
made into a minimal network. 

While reconstruction is not possible in general even for minimal networks 
(see Theorem 17.31 below) . the case of the torus teaches us that the map from 
conductances to the reponse matrix has potentially interesting preimage. 

A network ^ on a surface S is minimal if, when lift to the universal 
cover of S, the network is minimal, that is, the lifts of the medial strands do 
not self-intersect and no two lifts intersect more than once. 

As an example, see Figure [20] below. 

Theorem 5.1. Every network is topologically equivalent to a minimal net- 
work. 

Proof. Any surface with a boundary (or closed surface of genus > 2) has a 
hyperbolic metric (that is, a Riemannian metric of constant curvature —1) 
since its universal cover is the Poincare disk. 

Using this metric we can isotope the strands, fixing the endpoints, us- 
ing a curve-shortening fiow until all strands become hyperbolic geodesies. 
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Through this isotopy each singularity encountered either reduces the num- 
ber of crossings or is a triple crossing, which corresponds to a Y-Delta move 
of the network. In the end strands which touch the boundary will have the 
property that, when lifted to S, they meet other strands at most once. More 
generally, strands in different homotopy classes meet each other at most once 
in the universal cover. 

If two or more intersecting strands with no endpoints (that is, which are 
closed loops on the surface) have the same homotopy class then we must 
be more careful since their geodesic representatives will be identical. The 
situation in the universal cover will resemble that of Figure [131 l^ft, that 
is, we have a packet of strands which lie near the same geodesic, and other 
strands cross this packet transversely and exactly once each. 

It remains to prove that in this situation we can replace the self-intersecting 
packet with a "combed" version of itself, where all strands are parallel and 
non-intersecting (they will still intersect the transversal strands). For this it 
suffices to consider the situation on an annulus in which the packet strands 
wind around the annulus forming closed loops, and a certain number of other 
strands cross this packet transversely. On the universal cover of the annulus, 
which is a bi-infinite strip, the situation is as on the left in Figure [131 

We consider, then, a periodic graph on a strip Mx [0, 1] with conductances 
invariant under x H- x + 1. Suppose we now take n large and truncate the 
strip, by removing everything far to the left of the origin, say left of — n, 
leaving free boundary conditions there, and removing everything to the right 
of +ra, and similarly leaving free boundary conditions there. What remains 
is a planar graph on a rectangle [— n, ra] x [0, 1], and for this planar graph 
we can "comb out" the packet, as in Figure [T3l from the left to the right, so 
that the strands of the packet become parallel paths, nonintersecting except 
just left of the right endpoint, where they may cross. 

We claim that on the resulting combed finite planar graph, the conduc- 
tances are nearly periodic, that is, a periodic function plus an error tending 
to zero as the distance to ±n increases. To see this we use the interpretation 
of the response matrix entries as ratios of weighted sums of spanning forests. 
In [23] it is proved that for fixed i,j the truncation only changes Ljj by 
amounts tending to zero as n — )■ oo. It is also shown that Ljj — )■ as the 
distance from i to j increases. 

Take node i on the top and j on the bottom boundary of the strip. The 
numerator of the entry Lij is the weighted sum of spanning trees wired at 
the boundary and with a component connecting i and j (see Theorem | 
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We claim that for a typical tree the path from i to j does not get near ±n 
with high probability, that is, this path does not get far from i or j. If the 
path extended far beyond j, then (letting j + k denote the translate of j by 
k) for large k Lij+k would be of the same order as Lij, since a local change 
of the long path from i to j would result in a path from i to j + k. This 
contradicts the fact that Lij^^ tends to zero for large k. 

Thus Lij is a function depending less and less (as n — t- oo) on the con- 
ductances far from i and j. This proves convergence of the conductances 
as n — >■ oo, and the translation invariance of Lij implies the translation 
invariance of the limiting conductances. 

D 




Figure 13: Combing a river with three strands. 



5.3 Cycle-rooted spanning forests 

Where as the standard Laplacian determinant counts spanning trees, in the 
case of a one- or two-dimensional vector bundle with connection the bundle 
Laplacian determinant counts objects called cycle-rooted spanning forests. A 
cycle-rooted spanning forest (CRSF) is a subset S of the edges of ^, with 
the property that each connected component has exactly as many vertices as 
edges (and so has a unique cycle). See Figure | 



Theorem 5.2 ([llj). For a line bundle on a finite graph, 
detA= ^ JJce Yl (2-w-l/w) 

CRSFsT eeT cycles of T 

where the sum is over all CRSFs T , the first product is over the edges of 
T , and the second product is over the cycles of T, where w, 1/w are the 
monodromies of the two orientations of the cycle. 



U 



Figure 14: A CRSF on a square grid. 



Notice that for a line bundle the monodromies depend only on the cycle 
not on the starting point of the cycle. 

Theorem 5.3 ([14J). For a £?-hundle on a finite graph with SL2{'C) connec- 
tion, we have 

qdetA= 5^ n^^ n (2-TrH) 

CRSFsT eeT cycles of T 

where the sum is over all CRSFs T , the first product is over the edges ofT, 
and the second product is over the cycles of T, where w is the monodromy of 
the cycle (starting from some vertex, and in some arbitrary orientation). 

Notice that the trace of the monodromy is independent of starting point 
(since conjugation does not change the trace) and orientation, since for a 
matrix M G SL2(C), we have Tr(M) = Tr(M-i). 

Here the function qdet is the quaternion determinant of the self- 
dual matrix A; this requires some explanation. A matrix M with entries 
in GL2(C) is said to be self-dual if Mij 



Mji, where 




d 
— c 
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Since Ajj = —CijCpij where 0jj G SL2(C), we have Aji = —Cijcj)^^^ = Aij and 
so A is self-dual. 



We define 



qdet(M) = Y, sgn(^) J] I' 



-tYlw) 
adSn cycles 



where the sum is over the symmetric group, each permutation a is written 
as a product of disjoint cycles, and ti{w) is the trace of the product of the 
matrix entries in that cycle. If we group together terms above with the same 
cycles — up to the order of traversal of each cycle — then the contribution from 
each of these terms is identical: reversing the orientation of a cycle does not 
change its trace. So we can write 

c 

qdet(M)= Y. {-^r^^'Xl^Ay^i), 

cycle decomps i=l 

where the sum is over cycle decompositions of the indices (not taking into 
account the orientation of the cycles), c is the number of cycles, and Wi is 
the monodromy (in one direction or the other) of each cycle. Here tr is equal 
to the trace for cycles of length at least 3; cycles of length 1 or 2 are their 
own reversals so we define ti{w) = |tr(w7) for these cycles. 

As an example, let A = aL C = cl and B = \ , , | . Then 

Qdet ( ^ ^ j = tr{A)tr{C) - ti{BB) = ac - {bih - &2&3). 

Note that if M is a self-dual n x n matrix then ZM, considered as a 
2n X 2n matrix is antisymmetric, where Z is the matrix with diagonal blocks 
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. and zeros elsewhere. The following theorem allows us to compute 
Q-determinants explicitly. 

Theorem 5.4 ([9J). Let M be an n x n self-dual matrix with entries in 
(7-^2 (C) and M' the associated 2n x 2n matrix, obtained by replacing each 
entry with the 2x2 block of its entries. Then qdet(M) = Pf(ZM'), the 
Pfaffian of the antisymmetric matrix ZM' . 

Note that the matrix A' (obtained by replacing entries in A with the 2x2 
block of their coefficients) is just the matrix A acting on the total space of 
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the bundle Wg. So up to a sign we can write 

qdetA = (detA)i/2_ 

5.3.1 Unitary connections and measures 

In the case of a hne bundle, if Hcycies 2 — w — l/w > for every CRSF, we 
can define, following Theorem I5.2[ a probability measure on CRSFs in which 
a CRSF has probability proportional to the product of its edge weights times 

the factor Hcycies 2 - w - l/w. 

For example in the case that \(j)e\ = 1 for every edge, that is, (f) is a. 
unitary connection, then A is a Hermitian, positive semi-definite matrix 
(and positive definite for generic unitary connections). In this case |w| = 1 
for all monodromies and so 2 — w — l/w > 0, and equal to zero only for cycles 
with trivial monodromy. 

Similarly in the case of a two-dimensional bundle if 0e ^ SU2 for all edges 
then A is unitary and 2 — Ttw > for all monodromies. 

In both cases we call /x = /i$ the corresponding probability measure. 

Examples of natural settings of these measures are when the graphs are 
embedded on surfaces with geometric structures, such as a Riemannian met- 
ric (in which case the Levi-Civita connection defines parallel transport of 
vectors in the tangent bundle across an edge and provides a ?7(l)-connection) 
or an SU2 structure on the surface and parallel transport restricted to the 
graph gives a SU2 connection. 

In the case of a flat connection, the monodromy of a contractible cycle is 
trivial, so that the probability measure is supported on CRSFs whose cycles 
are all topologically nontrivial. We call such CRSFs essential CRSFs. 

5.4 Cycle-rooted groves 

For the bundle Laplacian, the natural objects replacing groves are cycle- 
rooted groves. Given a graph Q embedded on a surface with nodes A/", a 
cycle-rooted grove (CRG) is a collection of edges with the property that 
each component is either 

1. A tree containing one or more nodes 

2. a cycle-rooted tree containing no nodes. 
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Theorems 15.21 and 15.31 have an extension to the case of a graph with 
boundary B G V; the results are that the bundle Laplacian determinants 
(for the Laplacian with boundary) are weighted sums of cycle rooted groves; 
the weights of components with cycles are as before, and the weights of tree 
components are just the product of their edge weights (there is no "mon- 
odromy" contribution for tree components). See {21j . 

6 Annular networks 

6.1 Laplacian determinant 

Let ^ be a network with fiat C*-connection on an annulus, and with no nodes. 
Let z be the monodromy around a loop generating the homotopy group of 
the annulus. Then any simple closed loop on Q will have monodromy 1 if 
contractible, z or z~^ if noncontractible. 

As a consequence by Theorem 15.21 the Laplacian determinant is 

oo ^ 

detA = Y,Ck{2-z~-)\ (16) 

fc=i 

where Ck is the weighted sum of CRSFs having k components. 

Thus P{z) = det A is a Laurent polynomial in z from which one can 
extract the information about the number of loops in a random sample of a 
CRSF on Q. For example substituting x = 2—z—l/z we have that P{x)/P{l) 
is the probability generating function for the number of components. 

Independently of the conductances, P{z) contains topological information 
about the graph Q, for example the highest power of z is the maximum 
number of disjoint loops which can be drawn on Q, each winding around the 
annulus. This is because any such set of loops can be completed to a grove 
by adding some edges. 

There is a simple characterization of the roots of P. By (I16p the roots 
come in reciprocal pairs one of which is a double root at z = 1. 

Theorem 6.1. P is reciprocal and the roots of P are real, positive, and 
distinct except for a double root at z = 1. 

Proof. Our first goal is to minimize the graph so that it is a string of loops 
as illustrated in Figure [HI each loop winding around the annulus. 
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As in the proof of Theorem 15. II above, we use the curve-shortening flow on 
the medial strands until we have only a packet of strands winding around the 
core of the annulus as in Figure [T3l left (but without the vertical strands). In 
the resulting network, choose vertices fi,f2, one adjacent to each boundary 
and temporarily consider these to be nodes. There is one strand wrapping 
outside Vi and one wrapping outside V2- Break these strands so they have 
two stubs near each of f i, ^2. Now continue to minimize the medial strands 
as in the proof of Theorem 15.11 In the minimal network, the stubs from Vi 
must cross to the stubs of V2 (rather than to each other), since otherwise the 
network would be disconnected. We can push all intersections of the finite 
strands (that is, those connecting stubs) to one side of the intersections of the 
stub/river strands, as described in Section 16.4.11 and Figure 1201 below. The 
network and medial graph now resemble that in Figure [161 The graph can 
now easily be converted into a string of loops using the sequence of operations 
illustrated in Figure [T71 




Figure 15: Every graph on an annulus can be minimized to a string of loops. 

Let ai, . . . ,an-i be the conductances on the edges of the string and 
bi, . . . ,bn the conductances on the loops. Then the Laplacian is a tridiagonal 
matrix (using X = 2 — z — 1/z) 

/«! + biX — Oi ' 

—ai ai + 02 + b2X —02 

A= -02 

O'n-2 + O.n-1 + ^n-lX —CLn-l 

\ -a„-i an-i + bnX.J 

We pre- and post-multiply by the diagonal matrix U whose diagonal entries 
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Figure 16: Intermediate stage in minimization of network on an annulus. 
(Medial strands dashed/red). 



are h, to get 



/t + ^ 



-Ql 



UMJ 



/b\b2 



\ 



-a.\ 



/bite 





t+t+^ 



V 



-a2 





Vtete 




an-2 1 Qn-l 
bn-1 bn-1 




-Ctn-l 



+ x 



yC 



•^bn-lbn 



Thus we see that for each root 2; of P, — 2 + 2; + - is an eigenvahie of 



/ 



M 



Bl. 

bi 

-ai 

Vbib2 


-ai 

Vbite 
0.1 _|_ 02 

b2 ~'~ b2 




— 122 












an-2 1 Qn-l 

bn-1 b„-i 

-In-l 



\ 



yjbn-lbn 
O.n-1 



\Jbn-\bn 



b„ 



J 



This matrix M is symmetric and so has real eigenvalues, which are positive 
by flT6|) . Symmetry also implies that M has orthogonal eigenvectors. Thus if 
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Figure 17: Passing from the intermediate stage to the string of loops: (1) 
From the network in Figure [161 use Delta-Y moves on the upper row of 
triangles. (l)-(2) Flatten the degree-2 vertices on the upper row and use 
Y-Deltas on the second row of vertices. (2)-(3) Use Delta- Ys on the row of 
triangles just above the squares. (3)-(4) Use Delta- Ys on the upper row of 
triangles. (4)-(5) Y-Delta. (5)-(6) Flatten and then Delta-Y. (6)-(7) Flatten. 



51 



—X is a multiple eigenvalue of M it has an eigenspace V^^x of dimension > 2. 
This means there is a nonzero eigenvector of M with eigenvalue —X whose 
first coordinate is 0. But the coordinates of an eigenvector (xi,a;2, . . . ,Xn) 
satisfy a length-three linear recurrence: 

ai fli 

AXi = —Xi -^===X2 

Ol V&1&2 

and for i > 1 

•. '^i— 1 /'^j— 1 <^ i \ i 

AXi — / Xi-i + [— h -J-jXi / , ^i+l- 

VOi-lbi Oi bi ^Jhihi+l 

Starting from xi = this implies (using that all Oj ^ 0) that Xi = for 
all i. This is a contradiction. We conclude that the eigenvalues of M, and 
therefore the roots of P, are distinct. D 

6.2 Cylinder example 

Let us compute, for a rectangular cylinder, det A for the flat line bundle with 
monodromy z. This will allow us to compute the corresponding distribution 
of cycles in a uniform random essential CRSF. Let Hra,n be the square grid 
network on a cylinder, obtained from the n x m square grid by adding edges 
from {n,i) to (0,i) for each i G [l,""!]. Hm,n is a "product" of an m-vertex 
linear network Qm and a circular network Z„ of length n. We put a flat 
line bundle structure on Hm,n by putting parallel transport z on the edges 
{n,i){0,i) and 1 on all other edges. 

The eigenvalues of the Laplacian Ag^ on the linear network are 2—2 cos — 
for k = 1,2, ... ,m (the corresponding eigenvectors are fk{x) = cos ^ '^^^ ' > 
for X = 0, 1, 2, ... ,m — 1). 

The eigenvalues of the line-bundle Laplacian Az„ on Z„ with monodromy 
z are 2 — ( — 1/( where ( is an nth root of z. (The corresponding eigenvectors 
are gk{x) = (^.) The eigenvectors on Hm,n are the products fkQi for {k,i) G 
[l,m] X [0,n — 1]. 

We then have 

det A$ = TT TT 4 - 2 cos — - C - -. 
J.J. J.J. ^ /• 
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Use the identity 

where Ch„ is defined by Ch„(a + -) = «" + a~" (a variant of the Chebychev 
polynomial). 
This leads to 



nn"TT 1 

(Ch„(4-2cos— )-2--) 



k=l 

m -, 

nnTT 
ti; + Ch„(4-2cos — ) -2 
m 

k=l 

where w = 2 — z ~ ^. If we let Q(w) be this polynomial then Q(w)/Q(l) 
is the probability generating function of the number of cycles in a uniform 
random essential CRSF. 

It is interesting to consider what happens to this distribution for a large 
annulus, when m, n — t- oo with m/n converging to a fixed quantity r. For 
large n, Ch„(4— 2 cos — ) is large unless k is near 0. Thus only values of k near 
affect the limiting distribution. We have 4-2cos^ = 2 + ^ + 0(-)^ = 
Uk + l/ofc where a^ = 1 + — + 0(— )^. Thus in the limit m,n — )■ oo with 
m/n — )■ r we have Ch„(4 — 2 cos — ) = 2 cosh — + o(l). The limit probability 
generating function for the number of cycles is then 



Q{w)/Q{l) = w\[ 



fc= 



w + 2 cosh ^ - 2 

T_ 

2 cosh ^ - 1 



See Figure 



6.3 The response matrix 



Lam and Pylyavskyy [24] studied the response matrices of networks on annuli. 
Their point of view was to consider the lift Q of the network Q to the universal 
cover of the annulus. There one can define the response matrix as the limit 
of response matrices on larger and larger portions of the graph; it is not 
hard to show that this limit exists. Since the universal cover is planar one 
can recover some of the results from the planar case like nonnegativity of 
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Figure 18: The probabilities of 1, 2, . . . , 6 cycles in a uniform random essential 
CRSF on a large annulus, as a function of the modulus r. 

the non-interlaced minors (non-interlacedness makes sense since the lift Q is 
planar) . 

We will take a different approach here which uses the bundle Laplacian 
for a flat connection. Let ^ be a network drawn on an annulus S, with nodes 
M which are a subset of the vertices adjacent to the two boundaries of S. Let 
$ be a flat C*-connection with monodromy z around a generator of 7ri(S). 
Let A be the associated bundle Laplacian. Let L be the response matrix, 
defined as before as 

L = -A + BC-^B* 

where A = I ^^ I • Now entries in L are rational functions of z (with 



^B* C^ 
coefficients which depend on the conductances). 

What can be said about EIT on annular networks? 

For circular planar networks it was very useful to have combinatorial 
interpretations (in terms of groves) of the entries and minors of the response 
matrix. We have similar interpretations in the present case. The following 
theorem of [2T] holds for any network with line bundle (not just networks on 
an annulus). Compare Theorem 14.31 above. 

Theorem 6.2 ([2I])- Let Q,R, S,T be a partition of N and \R\ = \S\. Then 
det L^^ is the ratio of two terms: the denominator is the weighted sum of 
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CRSFs; the numerator is a signed weighted sum of cycle-rooted groves ofQx, 
the graph Q in which all nodes in T are considered internal, the nodes in 
Q are in singleton parts, and in which nodes in R are paired with nodes in 
S, with the sign being the sign of the pairing permutation, times the parallel 
transports from S to R: 

p.rsp(i)| |^p(fe)|„ I |„ 1 

det Lj^ = y {-ly ^'' '^;'7 '^^ •••'^^^ 

^'^ Z^ ^ ^ Z\l\2\---\n] 

permutations p l i i i j 

Here the numerator ^[ri'^'l ■ ■ ■ Irfc'*"'!!?!! ■ ■ ■ ji?^] is the weighted sum of 
CRG's in which there are tree components connecting r^ to Sp(j) for each 
i] the weight is the product of the edge conductances, times the product over 
all cycles 7 of 2 — w{'~f) — l/w{'~f) (where 1^(7) is the monodromy of the cycle 
7) and times the product of the parallel transports from the Tj to the Sp(i) 
along the edges of the trees. We use the notation Z as opposed to Z to 
remind us that we are dealing with the bundle Laplacian. 

On an annulus with fiat connection, it is convenient to have the connection 
supported on a "zipper" , that is, the set of faces crossing a shortest path in 
the dual graph from one boundary component to the other, as in Figure [191 
In such a case the monodromy along a path from a node to another node is 
l,z or z~^ if the nodes are on the same boundary; it is a power of z if the 
nodes are on different boundaries, and it is z^^ for a topologically nontrivial 
cycle. 

If particular if i and j are nodes on opposite boundaries then 

where Ck is the weighted sum of CRGs with a component tree connecting i 
to j, no cycles (a cycle is precluded by the existence of a path from z to j), 
and where k is the signed number of times the path from i to j crosses the 
zipper. 

Thus LijZ^-p^c is a Laurent polynomial of z with nonnegative coefficients. 

Here is another example. Suppose A = {ai, 02} are nodes on one bound- 
ary and B = {bi, 62} are nodes on the other. Then by the Theorem L^Z^^^ 
contains terms in which ai is paired with 61 (and 02 with 62) and, with an 
opposite sign terms in which ai is paired with 62 and 02 with bi. However 
these two types of terms differ by an power of z. Suppose for example that 
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Figure 19: The parallel transport on an annulus can be supported on a zipper. 

the zipper is as shown in Figure [191 Then any term ^^^1^^ j^as an even power 
of z and ^^|^i has an odd power of z. Thus L^Z^-^c is a Laurent polynomial 
in z with coefficients of alternating sign. 

6.4 Minimality and reconstruction 

6.4.1 Minimal networks on the annulus: landscapes 

Let ^ be a minimal network on an annulus S with n\ nodes on one boundary 
and ^2 on the other. As discussed, let Q be the lift of Q to the strip, the 
universal cover of the annulus. The medial strands come in four types: 

1. they connect stubs on the lower boundary of the strip (the rocks); 

2. they connect stubs on the upper boundary of the strip (the clouds); 

3. they connect stubs across the boundaries (the trees); 

4. they form bi- infinite paths along the strip (the river). 

We can perform Y-Delta transformations so that the strands of the first, 
second and fourth type do not intersect (that is, no strand of the first type 
intersects a strand of the second or fourth type, although it may intersect 
other strands of the first type, etc.) As long as there is at least one tree 
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strand, then we can "untwist" the river strands so that they are all parallel 
and noncrossing. 

The picture is then like a landscape; the rocks at the bottom, interspersed 
with trees which grow up to mingle with the clouds, and the river runs 
through the trees but does not touch the clouds or rocks. See Figure |20] and 




Figure 20: A landscape: the general form of a minimal network on an annulus 
(we show only the medial graph, and on the universal cover of the annulus). 




Figure 21: Bow falls, Banff. 



6.4.2 Reconstruction 

The reconstruction of the conductances on a minimal annular network as a 
function of the response matrix is open in general. Lam and Pylyavskyy [21] 
worked out the reconstruction for the "grid" network of Figure [221 whose 
medial graph is a grid. 

The new feature of this reconstruction problem is that the solution is not 
unique. If the medial graph of Q has k topologically nontrivial cycles, then 
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Figure 22: The "grid" network on an annulus. 

there are generically k\ solutions to the reconstruction. The idea is that given 
a set of conductances on Q one can permute the parallel strands of the medial 
graph and get another solution. How does one permute two strands? This 
follows from the discussion of combing in Theorem 15.11 above. See Figure 
[T3l One introduces a pair of edges of opposite conductance c, — c; this has no 
effect on the current flow or response matrix. Now perform Y-Delta moves 
around the annulus to exchange the adjacent strands. There is a unique 
choice of c for which the added edge, when it comes back around, will again 
have conductance c; it then can be removed along with the parallel edge of 
conductance — c. 

In [21] it is conjectured that reconstruction is possible on all minimal 
annular networks. 

7 Periodic networks and networks on the torus 

Let ^ be a network on a torus. By this we mean a network Q, with no 
boundary, embedded on a torus in such a way that every complementary 
component is contractible. 

For concreteness we suppose the torus is T^ = M^/Z^ (but the network 
can be arbitrary) . The lifted network Q on the universal cover M^ of T^ is a 
periodic planar network. 

This case differs from the previous cases of circular planar and annular 
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Z2 



Z2 



Z\ 



Z\ 



Figure 23: The 2x2 grid on a torus; in this example all edges have conduc- 
tance 1 except for the one indicated which has conductance 3. 

networks in that the underlying surface has no boundary, so we do not in- 
troduce any nodes. Remarkably we can still have a complete theory about 
reducibility, minimal networks, and the reconstruction problem. 



7.1 The spectral curve of A 

Let $ be a fiat line bundle on Q with monodromy zi, Z2 around the standard 
generators of the homotopy group of T^. Let A = A$ the associated Lapla- 
cian. See for example Figure [231 In this example the associated Laplacian 
determinant P{zi,Z2) = det A is (with the vertices in the order indicated) 



Pizi,Z2) 



det 



/ 



v 



6 

-3 — zi 

-1-Z2 





-3 - z^^ 
6 


-1-Z2 



2 3 42:1 42:2 
3z-,+^T- Azi Z2 

Z2 Zi 



4 



4 - .,- 



4 
-l~z^ 

76 







2n 



-1 

-1-2 

4 



\ 



-1 
1 



52 



zl 



Z1Z2 



-7Qzi + z^ + ^-52z2 + 264. 



Zl 



Z2 



(17) 

In general P is a (Laurent) polynomial with coefficients which are polyno- 
mials in the conductances. It is called the characteristic polynomial of A. 
The curve {(21,22) G C^ | P(zi,Z2) = 0} is the spectral curve of A. Note 
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that P is symmetric: P{zi,Z2) = P{z^'^, Z2^). This is because we can write 
A = d*Cd where d* is the Hermitian transpose of d when \zi\ = \z2\ = 1. 

The Newton polygon of the network Q on the torus is the Newton 
polygon of P, that is, the convex hull of 

{(i, j) e Z^ I zl'zi^ is a coefficient of P}. 

It is a centrally symmetric polygon centered at (0,0). 
By Theorem 15.21 we can write 



CRSFs 



P{ZI,Z2) = Yl \llcA{'^-W~l/w) 

a,s(2-44-^rV)', (18) 

(r,s)^(0,0) 

where in the first line we used the fact that in a CRSF on a torus all cycles 
have the same homology class, and in the second line we sum over all homol- 
ogy classes (r, s) G Z^ (we need only sum over one of each pair (r, s), (— r, —s)) 
where {r,s) = {ik,jk) with i,j being relatively prime. The coefficient Cr^s is 
the weighted sum of CRSFs with k cycles of homology class {i,j), that is (if 
we orient the cycles in the same direction) of total homology class (r, s). 

This shows that the z[z2 coefficient of P{zi, z^) is only due to the CRSFs 
with cycles in direction r / s. Thus the Newton polygon N tells us which 
homology classes of CRSFs are possible for Q\ equivalently, the maximal 
number of vertex-disjoint cycles in Q which have homology class (i, j) is 
determined by the integer point in A^ which is farthest from the origin in 
direction (i,j). 

In the above example f lT7|) . P can be written in "Newton polygon form" 

1 
-4 -52 -4 
P = 3 -76 264 -76 3. 
-4 -52 -4 

1 

From this, one can write P in a unique way in form fllSp : 

P{zi,Z2) = 3{2-zi--y + 6A{2-zi ) + {2-Z2 )'+48(2-Z2 ) + 

Zi Zi Z2 Z2 

+ 4(2-^1^2 )+4(2----). (19) 

Z1Z2 Z2 Zi 
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Thus there are, for example, exactly 48 CRSFs with a single cycle which 
has homology class (0, ±1). 

7.2 UST on the plane 

Let us say briefly a few words about the UST on the plane. 

Let Q be the periodic planar graph which is the lift of ^ to M^. Pemantle 
showed that there is a UST measure on spanning trees of Q which is the 



limit of the UST measure on larger and larger tori Qn = Q jnl? . 
The Free Energy of this measure is the limit 

F = lim —\ogZn, 

n— >co n 

where Zn is the partition sum for spanning trees on Qn- 
There is a nice formula for F 



{^T^^Y J\z^\=lJ\z2\=l Zl ^2 



'\Z1\ = 1J\Z2\ 

which arises from the exact formula 

Zr,{z,,Z2)= n n ^(^l'^2) 

(obtained by Fourier analysis of Qn) by taking logs and replacing the Riemann 
sum with an integral. 

One can analyze the limiting UST measure using the fact that edge pro- 
cess is determinantal with kernel T given by the Transfer Current matrix: 
T(e, e') is the current across e' when one unit of current enters the graph at 
e+ and exits at e_. In terms of the Greens function we can write 

T(e, e') = G(e+, e'+) - G(e+, e'_) - G(e_, e'+) + G(e_, e'_). 

The Green's function G{v,v') on Q has a formula involving the Fourier 
coefficients of 1/ P{zi,Z2): more precisely, an integral over the unit torus 
ki| = \z2\ = 1 of an expression z^Z2Q{zi,Z2)/P{zi,Z2) where (x,y) is the 
translation from v to v' and Q is a polynomial depending on where f , v' sit 
in their respective fundamental domains. In particular the transfer current 
has an exact asymptotic expression. 
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7.3 Medial strands and minimality 

Recall that two networks on T^ are topologically equivalent if one can be 
obtained from the other (disregarding conductances) by electrical transfor- 
mations. We say that they are electrically equivalent if their characteristic 
polynomials are the same (compare this definition with the one for networks 
with boundary in section l4.4.2p . 

Theorem 7.1. Two networks on a torus are topologically equivalent if and 
only if their Newton polygons are equal. 

Proof. It is not hard to show that electrical transformations do not change 
N. The simplest reason is because there is a bijection on the CRSFs of the 
"before" and "after" networks, and the set of CRSFs determines N. 

So it suffices to show that N determines a minimal network up to electrical 
equivalence. Suppose N has 2n integer boundary points vq, . . . , V2n-i in cclw 
order. Let Cj be the edge joining Vi to Vi+i, it is antiparallel with en+t 
because P is centrally symmetric. The proof now follows from the following 
lemma. D 

Lemma 7.2. The strands of a minimal network are in bijection with, and 
have homology classes equal to, the n edge pairs {ej,ej_|_„}. 

Proof. Since the network is minimal we can isotope the strands on the torus 
so that they are geodesies for a Euclidean structure. Then we orient the 
strands so that their y coordinate is increasing, or if they have no y compo- 
nent, the X component is increasing. 

Let Si, . . . , Sk be the strands, and let e, = (xj, yi) G Z^ be the homology 
class of Si. The length 2fcj of Si (that is, the number of edges of the corre- 
sponding zig-zag path of Q) is then determined by 2fcj = ^ • |ej A ej\ since 
exactly two strands cross at each edge of Q. 

Every strand gives a constraint on the homology class of a CRSF: for any 
CRSF with homology class cu we must have 

\u ACil < ki, (20) 

since |a; A ej| is the intersection number of the CRSF with the strand Si. 

Suppse the Si are indexed in order of decreasing slope Xi/yi. Let A^^ be 
the convex polygon with side vectors ei, . . . , e^, — ei, . . . , —Ck in cclw order. 
It has vertices vi,V2, . . . ,Vk, — f i, — f2, • • • , —Vk where 

.4(e......e,_.-e,-...-.). 
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We claim that A^^ is the subset of w G M^ which satisfy the constraints (120]) . 
To see this, note that if fi,fi+i are consecutive vertices of Ns and Cj is the 
vector between them, then 

2viAei = ei A Cj H h ei_i A e^ - Ci+i A e^ e^ A Cj 



E 



Cj A Cj I 



< 2ki 

and similarly for fj+i. Thus fj,fi+i satisfy the maximal possible constraint. 

Thus A^ C Ns- To see that A^ = Ns, is suffices to construct a CRSF with 

homology class equal to each of the vertices Vi of A^^. To construct a set of 

cycles with homology class Vi = |(— 61—62 ■ — 6^), for example, put weight 

— 1/2 on each edge of the zig-zag path of Si (in the orientation defined above). 
Each edge of Q will then have weight or ±1; since we assumed that the 
strands were drawn geodesically, at each vertex there is exactly one incoming 
and one outgoing edge. The oriented edges with weight 1 form a set of vertex 
disjoint cycles of the homology class u. D 

7.4 Harnack curves and characterization of spectral 
curves 

The Laurent polynomial P = det A has the following properties: 

1. P has real coefficients. 

2. P is symmetric: P(zi,Z2) = P{z^^, Z2^). 

3. P(l,l) = 0. 

4. {P = 0} is a Harnack curve. 

What we call Harnack curves were studied by Harnack and called simple 
Harnack curves by Mikhalkin in [26]. The definition given by Mikhalkin in 
[26] involves a topological condition on the real locus; in [27] a more concise 
characterization was given: a curve is Harnack if it intersects each torus 
{|^| = Ci, |u7| = 62} in at most two points (which are necessarily complex 
conjugates if not real). In particular the map {z, w) i-)- (log \z\, log \w\) which 
maps the curve into M^ is generically 2 to 1. The image of P = under 
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Figure 24: The amoeba (shaded) of the curve of Figure [23] above. The 
boundary of the amoeba is the real locus of P = 0; a model of P = can 
be obtained by gluing two copies of the shaded region together along their 
common boundary. 



this map is the amoeba of P. See Figure [23] for the amoeba of the example 

above. 

Remarkably, the four properties above characterize the characteristic poly- 
nomials of networks on a torus: 

Theorem 7.3. [151 HM -^'^'^ ^^^2/ polynomial P satisfying the above properties 
there is a network Q whose characteristic polynomial is P. The set of minimal 
networks with polynomial P , modulo electrical transformations, is a real torus 
of dimension g/2, where g is the geometric genus (the number of real ovals 
of the amoeba of P). 

8 Other surfaces 

For surfaces other than the annulus and torus, the fundamental group will 
be nonabelian so it is appropriate to use an SL2(C)-connection on Q rather 
than a C*-connection. 

Although we understand very little about this situation, there is one 
positive result. Considering the Laplacian determinant as a function on the 
moduli space of fiat SL2(C)-connections, it can be written in a unique way 
as a sum over CRSFs according to the homotopy types of their cycles. Let 
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us explain this. We define a finite lamination on E to be an isotopy class 
of finite pairwise disjoint collection of simple closed curves, none bounding 
a disk. For example on an annulus there is a finite lamination for each 
nonnegative integer k, k being the number of its cycles. 

The cycles in a CRSF form a finite lamination. One can ask the question: 
for a given finite lamination C, what is the weighted sum of CRSFs having 
cycle set of type £? We have 

Theorem 8.1 ([III [ID]). P = det A can be written as 

c 

where the sum is over all finite laminations and Cc is the weighted sum 
(partition sum) of CRSFs of cycle type C. The coefficients Cc are functions 
of P only, and can be extracted via an appropriate integration of P over the 
representation variety Hom^iiiiT,) , SL2{'C)). 

As an example, take S to be a pair of pants, that is, a sphere minus three 
disks. Let ^ be a network embedded on E. Take a fiat SL2(C)-connection on 
Q with monodromy A,B,C around the three boundary holes, where ABC = 
1. 

Then 

P = det A = Y^ CM.fc(2 - TrA)*(2 - TiBy{2 - Tt{AB))''. 

i,j,k>0 

Given variables X,Y,Z one can choose matrices A,B in SL2(C) such that 
2 - TrA = X, 2 - Trfi = F, 2 - TiAB = Z. Then 

P = P{X, Y,Z)= Y, c,,j,kX'Y^Z\ 

i,j,k>0 



One can extract from this expression the coefficients c 



i,j,k- 
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